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Abstract. We consider self-similar potential flow for compressible gas with 
polytropic pressure law. Self-similar solutions arise as large-time asymptotes 
of general solutions, and as exact solutions of many important special cases like 
Mach reflection, multidimensional Riemann problems, or flow around corners. 

Self-similar potential flow is a quasilinear second-order PDE of mixed type 
which is hyperbolic at infinity (if the velocity is globally bounded). The type 
in each point is determined by the local pseudo-Mach number L, with L < 1 
resp. L > 1 corresponding to elliptic resp. hyperbolic regions. We prove an 
ellipticity principle: the interior of a parabolic-elliptic region of a sufficiently 
smooth solution must be elliptic; in fact L must be bounded above away from 
1 by a domain-dependent function. In particular there are no open parabolic 
regions. We also discuss the case of slip boundary conditions at straight solid 
walls. 



Mixed-type equation; self-similar flow; compressible potential flow. 

1. Introduction 

We consider the isentropic Euler equations of compressible gas dynamics in d 
space dimensions: 

(1) pt + V • (pv) = 

(2) (pv) t + V ■ (pv®v) + V(p(p)) =0. 

Hereafter, V denotes the gradient with respect either to the space coordinates 
) or to similarity coordinates t l x. Throughout this paper, su- 
perscripts represent vector components whereas subscripts represent partial deriva- 
tives, v = (v 1 , v 2 , ■ ■ ■ , v d ) is the velocity of the gas, p the density, and p = p(p) 
follows the polytropic pressure law 



(3) P(P) 



CqPo / P 



7 \Po 



for 7 7^ (here cq is the sound speed at density po). 

For smooth solutions, substituting into J2J yields the simpler form 

Id) v t + Vt? • v + V(tt(/o)) = 0. 

Here tt is defined as 

<4) ^-{A^r-'-^M. 7*1 

l c olog^. 7=1- 

The Euler equations possess highly unstable vortex sheets. To focus on the 
acoustic waves and shock waves, one often considers the irrotationality assumption 
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Vj — v\ (where i,j = 1, . . . , d), which leads to 

v = V s 4> 

for some scalar potential function 4>. For smooth flows, substituting this into f2*Tjl 
yields, for i = 1, . . . , d, 

Q = (f>it+ ■ + n(p)i 

Thus, for some constant A, 

Substituting this into Q yields a single second-order quasilinear hyperbolic equa- 
tion, the potential flow equation, for a scalar field <f>: using c 2 = p p and (7r -1 )' = -4-, 
we get 

d 

(5) <j> tt + 2V& • V0 + ^ 0i0j<fe - c 2 A0 = 0. 

Due to the invariance of first-order systems of conservation laws under the scaling 
x <— sx, i <— si for s > 0, there is an important class of solutions to Q and (J2U, 
the self-similar solutions, which depend only on the similarity coordinates 

Self-similar solutions arise as large-time asymptotes of unsteady solutions and as 
exact solutions to many important special cases like Riemann problems, Mach re- 
flection, flow around solid corners etc. 

For potential flow the appropriate self-similar ansatz is 

4>(t, x) = #(|) 

because v — = V|^. (0 turns into the self-similar potential flow equation: 

d 

(6) c 2 A^ - Y, (ifc - €)&i ~ f>y = 0. 
The convenient change of variables 

, 111 2 

yields an equivalent equation whose coefficients do not directly depend on £: 

d 

(7) c 2 A x - J2 XiXiXH = IVxI 2 - dc 2 

For polytropic pressure laws, the sound speed c has a particularly simple form. 
For 7=1 (isothermal flow), c = co; for 7^1, 

(8) c 2 = ( 7 -l)(A-x-^), 
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Clearly, adding a constant to tp resp. x does not matter in © resp. Q, so we may 
take A = in (jHJ for simplicity: 

(9) c 2 = (l- 7 )(x+^). 

Vy = Vi/> — £ is called pseudo-velocity. 

Remark 1. (0 inherits a number of symmetries from |T|. (|2"|): 

(1) it is invariant under translation (which corresponds to the transformations 
v <— w + wo, x^x — «ot (mf/i parameter vq £ R d J in (t, x) coordinates), 

(2) it is invariant under rotation, and 

(3) it is invariant under velocity scaling (i.e. if x * s a global solution, then 
x(0 ; = s 2 x( s £) * s a global solution as well, for any s > 0). 

For many interesting flow patterns, 0) is of mixed type. The type is determined 
by the (local) pseudo-Mach number 

L :=!**!, 

c 

with < L < 1 for elliptic (pseudo-subsonic), L = 1 for parabolic (pseudo-sonic), 
i > 1 for hyperbolic (pseudo-supersonic) regions. 

It is important to identify the elliptic and hyperbolic regions because their math- 
ematical and physical properties are quite different. The equations 0) and are 
highly nonlinear, so the regions are not known a priori from the flow data. In con- 
trast to steady potential flow which has many interesting solutions of a single type, 
any self-similar potential flow solution in the entire space with bounded velocity 
must be hyperbolic at infinity and most cases appear to have an elliptic region as 
well. 

Our main result, Theorem states that if the local pseudo-Mach number of a 
C 3 solution of 10 is < 1 on the smooth boundary of a bounded domain fi, then it 
is bounded above away from 1 in the interior by a barrier function which depends 
only on the domain and on an upper bound on the sound speed. Such an example 
is useful for a variety of reasons; in the construction of self-similar solutions it guar- 
antees that perturbing the smooth boundary of a bounded elliptic region in a way 
that keeps L < 1 on the boundary must 1 keep the solution elliptic inside. In ab- 
sence of such a result, a possibly infinite number of parabolic- hyperbolic "bubbles" 
could arise in the interior which would make analysis prohibitively complicated. 

A variety of maximum principles for fluid variables has appeared in the literature. 
A direct analogue is the ellipticity principle derived in |Zhe97| (see also [LZY98 
IZhe01| ) for the self-similar pressure-gradient equations. [$cr91 shows that, under 
certain conditions, pressure in steady inviscid incompressible flow has extrema only 
at rest points. 

In [GT831 Chapter 15], rather general maximum principles for the gradient of 
second-order quasilinear elliptic PDE are provided. They apply in particular to 
equations of the form 

a ll {Vu)d lj u + b l {Vu)d l u = 0, 

4t is expected that for boundary conditions that admit a unique solution, the solution is C 3 
(in fact analytic) in the interior and that continuous perturbations of the boundary and boundary 
data yield continuous changes of the solution in any C k norm. 
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i.e. in absence of Oth order terms and with coefficients depending on the gradient 
only. Among the corollaries are strong maximum principles for the Mach number 
(or equivalently density or velocity) in steady potential flow as well as for the 
pseudo-Mach number L in the isothermal case (7 = 1) for self-similar potential 
flow. However, the case 7 7^ 1 is not covered. 



2. The ellipticity principle for self-similar polytropic potential 

FLOWS 

Theorem 1 (Ellipticity principle). Let d > 2 and 7 > — 1. Let f2 C M. d open and 
bounded. There is a positive constant S so that, for any c > and any b € C 2 (f2) 
with \Vb\ < §, |V 2 6| < Jr and for any solution X S C 3 (n) of with L < 1, p > 
and c < c, either 

L 2 < 1 - 5 



Remark 2. At its heart the bound depends only on the domain, since the velocities 
c. v, £ can be rescaled (see Remark^) to allow c — 1. 

Proof. During the proof, all equations and inequalities are meant to hold in the 
maximum point only. The O notation is with respect to S J, and uniform in c. 

Assume that L 2 + b has a maximum in some interior point with L 2 > 1 — S. We 
fix S < 1 right away, so we may assume |V%| = Lc =/= 0. Since the equations are 
rotation-invariant, we may, without loss of generality, choose xi > and \j = 
for all j > 1. The first-order conditions for a maximum imply 

= ^ + 6), = 2X1X11 ° 2 + (7 " + * u) + h 

c 4 

-c6i - (7 - l)L 3 I-7 
" Xll= L(2 + (7-l)L 2 ) = 7TT + Q((5) - 

For j > 1, 

= (L 2 + b) J = 2XlXljC2 + ( 7' 1)x ^+b 2 



^ Xij = 0((5). 

From 0, 



(10) 



5^ Xij = (i 2 - l)(xn + l) + L 2 -o!=l-d + 0(6) 
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Moreover, using (//<?)" = (f'g 2 ~ fgg" + 2f(g') 2 - 2fgg')/g 3 , the second-order 
conditions yield 

0> (L 2 + 6)n 

_ 2(xiXm +Xn +E J >iXi 3 )c 4 



(7 - i)xic 2 (xn(i + xn) + Ej>i xij + xixm) 
+ c 6 

2( 7 -l) 2 x j(l + Xll ) 2 4( 7 -l) x 2 c 2 xii(l + Xii) , _ anm 

+ ^6 + ^6 + C °W 

= c- 1 L(2+( 7 -l)i 2 )xm 

+ c- 2 (xn + 2( 7 - 1)L 2 (1 + xii))(2xn + (7 - l)i 2 (l + Xn)) 
+ c-\2+{ 1 -l)L 2 )Y J xl+c- 2 0{5) 

xin < cL(2 + ( ^_ 1)£2) ((xn + 2(7 - i)£ 2 (i + xii)K2xn + (7 - i)£ 2 (i + Xn)) 

=0(1) =0(<5) 
+ (2 + ( 7 - 1)L 2 )]T x 2 , + 0(5)) = c- l O(S) 

=0(S) 

For J > 1, 

> (L 2 + % 

= 2(xiXijj +EtiX 2 j)c 4 (7 - l)XiC^06j + Eti X| + XiXijj) 

L(2 + (j-l)L 2 ) (2 + ( 7 -l)£ 2 )(l + 2(7-l)L 2 ) 2 
~ Xijj + - 2 Xij 

2 /„, i\r2 



2+ ( 7 -i)^ 2 ( 7 -l)L _ 2 
+ ^ Z^Xkj + ~2 Xn+c 0{S) 



c Xljj < -1 ( (l + 2(7-l)£ 2 )x 2 J +E^)- 2 + Xjj+ ° {S) 

=-l+0(«) "" =( 7 -l)/(7+l)+0(5) 



c^Xi« < ~^X 2 



3>l j>l 



(11) + 2_l( d _i ) + ow 

7 + 1 

Minimizing X^>i X 2 j with respect to (see (1101) ^) 

= l-d + 0(J); 



the minimum is 
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l-d + 0(S) , 

X22 = ■ ■ ■ = xdd = — d _i = - i+ °( s y> 



substituting this into (|llfl yields 

(12) "Exyj^^ + oW 

At this point it is necessary to take di of JJJ: 



XjXfcXijfc 



>0 Kc-WiS) J>1 =0 ■ 7 ^ fe =0 



+ ((! - 7)xi + (! - 7)(xi Xn +E Xl i) ~ 2x1 Xl1 ) Xn 



=i^+0(5) J>1 =0 =^+0(5) =0(«)X1 



=0(8)xi 

- 7)Xi + (1 - 7)(xi Xn + X> ~ 2 y X^Xifc)Xfcfc 

fe>1 = ^+0(5) 3>1 ^ =0 

2^(xiiXi +XiXij) Xij -2 ^ (Xij Xk + Xj Xik)Xjk 
^8) 1<j<k ^ 

: rf(7-l)xi + (2 + rf(7-l))(Xi Xn 

Xk Xlk) 



^+o(s) k>1 



7 i>i i>i 

Using (fTOfl we see 

(13) c^ Xl „>0(^) 

i>i 

Comparing (|12|l and 113(1 yields a contradiction (for d > 2) for <5 chosen sufficiently 
small. □ 

Corollary 1. For d > 1, Q does noi admit solutions ^eC 3 £/iai are parabolic on 
some open set. 

Remark 3. Corollary^ does not hold for d = 1. In £/ws case Q reduces to the 
second-order ODE 



(14) (c^-(x') 2 )x" = (x') 2 -c 2 , c 2 



j(l-7)(x+i(x') 2 ), 77^1 
1 const, 7=1 



Solutions of the first-order ODE c 2 = (x') 2 satisfies (|14|l as we/?, 6m< i/iey are 
parabolic everywhere. (Note that those solutions correspond to rarefaction waves.) 

Nevertheless, solutions to (|14fl are affine on any non-parabolic interval (this is 
obvious from (jBJ for d — 1), so a weaker version of ellipticity holds trivially. 
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3. Straight solid walls with slip condition 

Theorem ^ by itself is suitable for parabolic boundaries (boundary condition 
L = 1) or regularizations (L = 1 — e). Here wc extend it to another common 
boundary condition, the slip condition 

< 15 » s-° 

at a solid wall. We consider straight walls only because in self-similar flow curved 
walls are not physically interesting. 

Theorem 2. Let d > 2 and 7 > — 1. Let U be a small open ball centered on 
{0} xR^ 1 and set Q := f7n((0,oo) xl") and T := U<1({0} xl^ 1 ). There is 
a positive constant 5 so that, for any c > 0, any b G C 2 (fiUT) with H = 
|Vo| < 4, |V 2 o| < 4r in O, and any x G C7 3 (fi U T) f/iaf sofoes Q) in f2, satisfies 
(|15|l on r and /ias L < 1, p > and c < c in flUT, either 

L 2 < 1 - S 

or L 2 + b does not attain its maximum in ft U T. 

Proof. In this setting (II 51) means Xi = on T. Taking tangential derivatives in 
(fT3f> yields 

Xh = 0, Xiw=0 Vi,j G {2, ...,d} 

on F. This implies 

(c 2 )i = (1 ~7)( Xi + Xi X11 + X^jXij) = 
=0 =0 J=2 =0 

on r as well. Moreover take d\ of (JJJ: 

d 

c 2 Xm + (c 2 )iXn + J2(° 2xiu + ( c2 )^Xu) ~ xi Xiu - 2 Xi X?i 
=0 1=2 =0 =0 =0 =0 

d d 

~ 2 'Y^(Xi_X]Xii3 + (XiXj)i Xij ) - (XiiXj + Xr Xij + XiXjXiij) 
3=2 =0 =0 1,0=2 =0 =0 =0 

d 



2 Xi Xil + 2 2^ Xj Xij ~d(c)i on T 
=0 3=2 =0 =0 



xm = on r 

Combining all results we see that even reflection 

x(e,e,---,t d ) --=x(-e,e,---,c l ) 

yields a C 3 (U) extension of x that satisfies (0 in all of U. Moreover since 61 = 
on r, we can extend b in the same manner to C 2 (U). All conditions of Theorem 
are satisfied, so it implies the desired result. □ 

Remark 4. Any solid wall can be rotated and translated (see Remark^) to coincide 
with r in Theorem^ 
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